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ABSTRACT 



Context. We have observed the Crab Pulsar in the optical with S-Cam, an instrument based on Superconducting Tunneling Junctions (STJs) 
with yus time resolution. 

Aims. Our aim was to study the delay between the radio and optical pulse. 

Methods. The Crab Pulsar was observed three times over a time span of almost 7 years, on two different locations, using three different 
versions of the instrument, and using two different GPS units. 

Results. We consistently find that the optical peak leads the radio peak by 49+90, 254+170, and 291±100 fis. On assumption of a constant 
optical lead, the weighted-average value is ~170 /js, or when rejecting (based on a perhaps questionable radio ephemeris) the first measurement, 
273±100/is. 
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1. Introduction 

Precise timing of pulsar light curves throughout the electro- 
magnetic spectrum is a powerful tool to constrain theories of 
the spatial distribution of various emission regions. In recent 
years, it has become clear that the main and secondary pulses 
of the Crab Pulsar (PSR J0534+2200) are not aligned in time at 
different wavelengths. X-rays are leading the radio pulse by re- 



ported values of 344+40 ^s (Rots et al. 2004 RXTE data) and 



280+40 yus (Kuiper et al. 120031 INTEGRAL data) and y-rays 
are leading the radio pulse by 241+29 fis (Kuiper et al. 2003 
EGRET data. The uncertainty in this value does not include the 
EGRET absolute timing uncertainty of better than 100 /vs). At 
optical wavelengths, the observations present a less coherent 
picture. Sanwal (1999) has reported a time shift of 140 /is (op- 
tical leading the radio). The uncertainty in this value is 20 /is 
in the determination of the optical peak and 75 yus in the ra- 
dio ephemeris. Shearer et al. (2003 1 have reported a lead of 
100+20 //s for simultaneous optical and radio observations of 
giant radio pulses. Golden et al. (2000) have reported that the 
optical pulse trails the radio pulse by ~80+60 /vs. Romani et al. 
(2001 1 conclude that the radio and optical peaks are coincident 
to better than 30 jt/s, but their error excludes the uncertainty 
of the radio ephemeris (150 jt/s). The internal inconsistency of 



these results - if we assume that the optical-radio delay is con- 
stant - has prompted us to look into this matter in detail using 
recent observations in combination with earlier data. 



2. Observations 

Observations were obtained with S-Cam3 in November 2005 
on the ESA Optical Ground Station (OGS) telescope on 
Tenerife and with S-Cam2 in October 2000 and S-Caml in 
February 1999 (see Perryman et al. Il999b on the WHT tele- 
scope on La Palma (Tabled- All S-Cam instruments are based 
on Superconducting Tunneling Junctions (STJs), which regis- 
ter individual photons with high time resolution and moderate 
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spectral resolution (see e.g. Martin et al. 2003 1. Times are pro- 
vided by a GPS unit which yields an accuracy better than 1 /zs 
with respect to UTC in individual photon arrival times. The 
time resolution of the instrument was 5 yus for S-Caml and S- 
Cam2 and 1 yus for S-Cam3. 

3. Analysis and Results 

The time tags of all photons were converted from UTC 
to Barycentric Dynamical Time (TDB) at the solar system 
barycenter (a process we refer to as "barycentering"). This pro- 
cess involves conversion from UTC to TAI (by - for our obser- 
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Table 1. Log of observations. 



Instrument 


Observation period (YYYY-MM-DD) 


Exposure (s) 


Observatory 


GPS Instrument 


S-Caml 


1999-02-02 - 1999-02-04 


4800 


WHT La Palma (4.2 m) 


HP 58503A 


S-Cam2 


2000-10-01 -2000-10-02 


3317 


WHT La Palma (4.2 m) 


HP 58503A 


S-Cam3 


2005-11-12 


2616 


OGS Tenerife (1 m) 


Zyfer GPStarplus Model 565 



vations - adding 32 leap seconds) and subsequent conversion 
from TAI to TT by adding the constant offset of 32.184 s. TT 
is then converted TDB. (Although the IAU officially replaced 
TDB by TCB in 1991, they reckoned that, where continuity 
with previous work is desirable, TDB may be still be used.) 
Furthermore, a correction from the observatory to the Earth 
centre and a correction for the gravitational propagation de- 
lay are applied. We have used the JPL DE200 ephemeris (and 
not the more recent DE405 ephemeris) to perform the bary cen- 
tering, since we compare our data against observations of the 
Crab Pulsar in the radio, which are also barycentered using 
this ephemeris. WGS84 coordinates of the telescopes were ob- 
tained with a GPS unit. These geodetic coordinates were con- 
verted into geocentric coordinates to obtain the proper obser- 
vatory to Earth centre correction. 

We have folded our data using the P, P, and arrival times 
(at infinite frequency) from the Crab radio ephemeris (Lyne et 
al.|2005l. The folded profiles are displayed in Fig. [2 An im- 
portant parameter from this ephemeris is the uncertainy in ra- 
dio timing, which we denote with cr ra( ji . This uncertainty is 
composed of different elements which can not all be treated 
as a statistical error that reduces when averaged over multiple 
observations. The most important components are unmodeled 
delays and imperfect polarization calibration (~40 fis), uncer- 
tainties in the delay due to scattering (~20 fis), and variations 
in the daily time-of-arrival measurements with rms residuals of 
the order 20-50 yus. 

No default method exists for determining the position of the 
peak of the profile. We follow the method of Rots et al. ( 2004 1, 
which differs, however, from the choice of other authors (e.g. 
Kuiper et al. 2003 1. The folded pulse profile (subdivided into 
1000 phase bins) was fitted in a phase range of 0.01 centered 
around the peak using an iterative approach, allowing the phase 
range in use to evolve with the fitted peak position. The model 
used to determine the peak was a Lorentzian with normaliza- 
tion, width, and center as free parameters. This choice is some- 
what arbitrary and different fitting functions (e.g. a Gaussian 
or a parabola) return slightly different results. Variations, how- 
ever, are always much less than the uncertainty in the radio 
ephemeris, the final number to which we compare against. We 
have also studied the behaviour of the peak position obtained 
with different phase ranges. These experiments show a consis- 
tent behaviour: due to the asymmetric shape of the pulse pro- 
file, the phase of the peak shifts to a higher value as the phase 
range is decreased but flattens off for phase ranges smaller than 
0.01. We estimate that the uncertainty introduced by using a 
finite phase range of 0.01 (instead of an ideal, infinitely small 
phase range) is 10 fis. Our results are summarised in Tabled In 
all cases, the optical pulse leads the radio pulse. The time shift 
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Fig. 1. The three renormalized profiles: S-Caml in black, S- 
Cam2 in red, S-Cam3 in green (in the electronic version). In the 
top panel, the total profile is displayed. The bottom panel shows 
a zoom-in on the main peak. The statistical noise is largest 
for the S-Cam3 observation, due to the smaller telescope and 
shorter exposure time. Phase (or 1) is defined as the peak 
of the radio ephemeris; this phase is indicated by the vertical 
dashed line. 



found from S-Caml data is clearly the smallest, while the value 
from S-Cam2 data is affected by a substantially larger cr m dio- 

We have also determined the position of the secondary peak 
using the same algorithm as above but with twice the phase 
range (i.e. a width of 0.02). This increased range is required 
since the secondary peak is less sharp (resulting in almost no 
curvature in a 0.01 phase range) and has a lower signal to noise 
ratio. The results (see Table |2j are consistent with a constant 
phase difference between the primary and secondary peak. This 
confirms that the S-Caml pulse profile is shifted (compared to 
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the S-Cam2 and S-Cam3 data points): not only the main peak 
arrives later, but also the secondary, indicating a shift of the 
whole pulse profile, not a change in the main peak only. 

In order to compare the results of the three S-Cam cam- 
paigns, the three different pulse profiles were scaled to match 
each other. The level in the phase range 0.7-0.8 was calculated 
as well as the peak level determined from the mean of the three 
highest bins. The profile was renormalized such that the peak is 
at 1 and the base level at 0. This scaling assumes that there are 
no significant counts from the pulsar in the 0.7-0.8 phase range 
and that all other photons (sky background, nebula background, 
and noise photons) are distributed equally in phase. 

The normalised profiles look similar, albeit the S-Caml 
profile looks a bit different at the peak: the leading part of the 
peak is almost identical to the S-Cam3 profile (up to phase 
0.98), but then starts to move away and peaks later than the 
S-Cam2 and S-Cam3 profiles. We have investigated the pos- 
sibility that the P and P which have been used to fold the 
data are sub-optimum by performing a period search on our 
data. However, the solution hence found is fully consistent 
with the radio ephemeris, with the latter having a higher pre- 
cision. Although we do not have evidence for a large colour- 
dependence of the profile shape, small differences in the de- 
tailed shape of the pulse profile may be explained by the differ- 
ent spectral response and array size of the different generations 
of the S-Cam instrument, possibly in combination with differ- 
ent sky-background levels and atmospheric conditions under 
which the observations were acquired. 

We have also investigated, by exploiting the intrinsic spec- 
tral resolution of STJs, whether there is a phase difference in 
the peak between the "red" and the "blue". We have selected 
the ~30% photons with highest energy (~400-500 nm) and the 
~30% photons with lowest energy (~600-760 nm). Because of 
the limited spectral resolution of the S-Cam instrument, this ap- 
proach yields a cleaner colour separation than splitting the pho- 
tons in two adjacent bands, each containing 50% of all events. 
We then determine the delay between the photons in the red 
band with respect to the photons in the blue band by doing a 
cross-correlation of the whole main peak profile. We find that 
the red photons lag by -1±2 //s for S-Caml, 4+4 fis for S- 
Cam2, and 5+10 fis for S-Cam3. This does not represent a sig- 
nificant time shift and can be compared to the results of Sanwal 
dl999> who found a delay between the U and R band of 10±4 //s 
and Golden et al. ( 2000 ) who found that the peaks in B and V 
are coincident within 10 jus. 

4. Checks on correctness of times 

The following checks were made: 

GPS time: It was verified that both GPS units provide the 
same time to within 0.1 fis, the typical intrinsic GPS time jitter. 

GPS position: The correct reading of the S-Cam3 GPS was 
verified by the GPS present at the OGS telescope. This GPS 
presents the position in Earth-Centered, Earth-Fixed (ECEF) 
coordinates (X,Y,Z) = (5 390284, -1597 899, 3 007075) m. 
The reading of the S-Cam3 GPS (latitude 28° 18. m 0640, lon- 
gitude 16° 30. m 7133 and altitude 2447.02 m) was converted 
into ECEF coordinates which resulted in (X,Y,Z) = (5 390281, 



-1 597 891, 3 007083) m. These coordinates differ by 12 m, 
consistent with the separation between the two GPS antennas. 

Period-folding software: Two different approaches of fold- 
ing the data were followed and it was verified that they gave 
identical results. The first approach is to fold the data with the 
period, period derivative, and epoch as given by the Lyne et al. 
(2005 1 radio-timing results. The second approach is to deter- 
mine, by phase folding, the "current period" (i.e. at the time of 
the observation) and fold the data on that period without tak- 
ing the period derivative into account, which is not necessary 
for such short datasets. The phase zero (or epoch) close to the 
time of the observation should be independently calculated us- 
ing the radio ephemeris. While in theory both methods should 
give identical results, this is by no means guaranteed in prac- 
tice: since the time of observation can be 15 days away from 
the ephemeris epoch, the calculation of the phase using the pe- 
riod and period derivative requires high precision. It was veri- 
fied that this requirement on precision is met by the FTOOLS 
"efold" program. 

Barycentering software: Using the detailed examples pro- 
vided in Lyne et al. (2005 1, it was verified that our barycenter- 
ing code works correctly at the 5 fis level (dominated by cor- 
rections in TT to TDB: relativistic corrections from the geoid 
to the solar system barycenter). 

Crab coordinates: For the barycentering of the Crab, the 
coordinates used by Lyne et al. are: 05 h 34 m 31. s 97232, +22° 
00' 52."0690 (J2000). For each monthly ephemeris, this fixed 
position is used as an assumed position for the pulsar at the 
relevant epoch. The known proper motion of ~18 mas/yr is ig- 
nored, and the effect of this proper motion on the true position 
is incorporated into the solution of the period, period deriva- 
tive and arrival time (C.A. Jordan, private communication). Our 
barycentric correction assumes the same reference position as 
the relevant monthly ephemeris, so that consistent arrival times 
are obtained. The erroneous use of proper-motion corrected co- 
ordinates in the barycentering would have a significant effect: 
inclusion of a proper motion of 18 mas/yr would imply a time 
difference of ~17 //sec per year, or ~ 100 yt/s for the 2005 obser- 
vation. 

Instrumental effects: The observation of delays/shifts at the 
level of a few hundred jus and the realization that this time scale 
is comparable to the average time between the events for the 
mean count rate during the Crab observations, raises the pos- 
sibility that photons could have the time tag belonging to the 
previous photon. This wrong association could either be caused 
by software or by instrumental effects. Both possibilities were 
studied in detail and rejected: we are confident that our photons 
have correct time tags. 

Instrumental effects (S-Cam2): The S-Cam2 data suffered 
from a hardware problem which resulted in corruption of time 
tags of 75% of all events, randomly distributed in time. Since 
this problem is fully deterministic and understood, photons 
with corrupted time tags have been filtered out. The remain- 
ing 25% of photons, which have reliable time tags, have been 
used for this analysis. 

Instrumental delays (S-Caml and S-Cam2): In S-Caml and 
S-Cam2, incoming photons generate electronic signals from 
which the photon energy and arrival time can be extracted. The 
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Table 2. Radio to optical delays. The used parameters from the Crab monthly ephemeris are also displayed (converted into P 
and P). The number of digits displayed is higher than the uncertainty. The cr radio quoted is the uncertainty in the radio ephemeris 
as published in the Crab monthly ephemeris. The separation is defined as the difference (in phase) between the primary and 
secondary peak. 



Instrument 


Delay (jis) 




Separation 


P(s) 


P(lQ- ,3 s s- [ ) 


Epoch (MJD) 


S-Caml 


49±10 


80 


0.4056±0.0009 


0.03349352792448 


4.20554737 


51224.000000352002 


S-Cam2 


254±8 


160 


0.4052±0.0005 


0.03351561813489 


4.20539152 


51832.000000238553 


S-Cam3 


291±13 


90 


0.4065±0.0017 


0.03358309316492 


4.20593301 


53689.000000194479 



Q 100 

o 

04 o 



S-Cam2 



Sanwal 



S-Caml 



Golden et al. 



S-Cam3 



Romani et al. 



Observation Date 

Fig. 2. Our three results combined with literature results in 
the optical. For all results, the uncertainties were obtained by 
linear addition (i.e. not quadratic) of the uncertainty in the 
peak determination and the uncertainty in the radio ephemeris. 
Uncertainties on the peak determination alone are also plotted. 



The open symbol refers to Shearer et al. (2003 1, who observed 
giant radio pulses. 



arrival time of a photon is defined as the moment in time that 
the electronic signal rises above a certain threshold. The delay 
of the registered arrival time compared to the true arrival time 
is less than 10 /l/s, and corrections for this delay have not been 
made for S-Caml and S-Cam2 data. 

Instrumental delays (S-Cam3): In S-Cam3, each photon 
gives rise to a bi-polar signal in the detector electronics chain. 
The zero-crossing of this signal is used to time tag the photon. 
The delay of this moment compared to the true time the photon 
entered the detector has been calibrated to be 66.0+ 1 .0 //s (fully 
consistent with the value of 66.7 fis obtained from simulations 
of the instrument). All S-Cam3 time tags have been corrected 
for this (fixed) delay. 

5. Discussion and Conclusions 

We have observed the Crab Pulsar with three different gener- 
ations of the S-Cam instrument, on two locations, using two 
GPS units. We consistently find that the optical pulse is lead- 
ing the radio pulse. However, the amount by which the opti- 
cal is leading the radio differs from observation to observa- 
tion. When comparing the different results as plotted in Fig. 
[5] which are obtained with the Jodrell Bank ephemeris, one 
should take into account that cr radio contains a systematical 



component (of ~40 fis) affecting all these measurements in the 
same way.If we subtract the 40 //s, for the purpose of compar- 
ing the results, from cr^ao and add the remainder, in quadra- 
ture, to our measurement uncertainties, we obtain the following 
values for S-Caml, S-Cam2, and S-Cam3: 49+41, 254+120, 
and 291+50yus. 

From the last two observations (S-Cam2 and S-Cam3), we 
then determine an average lead of 273±65 //s. The uncertainty 
in this determination is dominated by uncertainties in the ra- 
dio ephemeris. The result in X-rays for the same two epochs is 
370+40 //s (Rots, private communication). 

From our data, we cannot rigorously exclude the possibility 
that the delay between the radio and optical peak evolves with 
time. The S-Caml result (49+90 jus; the uncertainty has been 
determined by a linear addition of the uncertainty in the peak 
determination and cr radio ) in particular is slightly deviant com- 
pared to the S-Cam2 and S-Cam3 data points (254+170 and 
291 + 100 //s). However, the radio observations at the times of 
the S-Caml observations have a high dispersion measure - a 
higher value has not been observed after February 1999. As a 
result, the radio-ephemeris uncertainty (80 /vs) might have been 
slightly underestimated: variations in the dispersion measure 
could introduce a somewhat higher uncertainty than the stan- 
dard 20 yus (Lyne et al. 1993). We note that Rots et al. (2004i 
have discarded their contemporaneous (X-ray) data point, be- 
cause the timing ephemeris is somewhat suspect (the GRO ver- 
sion has a high second derivative indicative of a questionable 
fit). In the following we will therefore only consider the S- 
Cam2 and S-Cam3 data points. We emphasize here that cor- 
rections for arrival times at infinite frequency (which are sub- 
stantial: ~0.6 s at 610 MHz) depend on the dispersion measure. 
The only (literature) data point which is clearly inconsistent 
with all observations is that of Golden et al. (|2000l. We have 
no explanation for this. 

Our value of the optical phase difference between the main 
and secondary pulse of 0.4054±0.0004 is not consistent with 
the X-ray value from Rots et al. (120041 of 0.4001 ±0.0002 peri- 
ods. This implies that the details of the pulse profile are differ- 
ent in X-rays and at optical wavelengths. 

Our time shift of 273+65 /is is somewhat smaller than, but 
consistent with, the time shift as obtained from X-ray measure- 
ments. A time shift of ~270 fis indicates that possibly (in a 
simple geometrical model ignoring relativistic effects) the op- 
tical radiation is formed ~90 km higher in the magnetosphere 
than the radio emission. Alternatively, the difference in phase 
of ~0.008 could be interpreted as an angle between the radio 
and optical beam of ~3°. 
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Ideally, a simultaneous radio-optical observation at high 
frequency should be performed. With an observation like this 
uncertainties resulting from corrections for interstellar scatter- 
ing are minimized, and the accuracy will effectively be limited 
by systematic effects (~40 /ys for Jodrell Bank). 
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